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We discuss the evolution of topological states and their orbital textures in the mixed valence 
compounds SmBe and YbBe within the framework of the generalized gradient approximation plus 
onsite Coulomb interaction (GGA+U) scheme for a wide range of values of U. In SmBe, the 
topological Kondo insulator (TKI) gap is found to be insensitive to the value of U, but in sharp 
contrast, Kondo physics in isostructural YbBe displays a surprising sensitivity to U. In particular, 
as U is increased in YbBe, the correlated TKI state in the weak-coupling regime transforms into a 
d-p-iype topological insulator phase with a band inversion between Yb-5d and B-2p orbitals in the 
intermediate coupling range, without closing the insulating energy gap throughout this process. Our 
theoretical predictions related to the TKI and non-TKI phases in SmBe and YbBe are in substantial 
accord with recent angle-resolved photoemission spectroscopy (ARPES) experiments. 


I. INTRODUCTION 

Topological insulators (TIs) present a fundamentally 
new state of quantum mattei^Ilf^, which harbors a 
non-trivial bulk electronic structure with gapless sur¬ 
face states on the boundaries. Experimental discov¬ 
eries followed soon after the first theoretical proposals 
in both two-dimensional (2 dP3EI 

and 3D material^ 
were made. Subsequently, many distinct classes of 
TIs have been proposed, including: topological Mott 
insulatorJ^ topological Anderson insulatorJ^ topo¬ 
logical crystalline in^lator^m and topological Kondo 
insulators (TKIs)P^^. A variety of broken symmetry 
TIs have also been proposed including topological an¬ 
tiferromagnetic insulator^^, topological hidden ordeJ^, 
and fully gapped topological super conductors.!^ Among 
these many possibilities, TKIs have attracted much re¬ 
cent attention because of their potential realization in 
and YbBeP^M^, some inconsistencies between 
the predicted and experimentally observed properties in 
these systems notwithstanding. 

Several lanthanide and actinide materials that contain 
localized /-electron have been predicted to be topologi¬ 
cally nontrivial. For example, actinide based 5/ materi¬ 
als are predicted to be TIs when onsite Coulomb interac¬ 
tion is included in the computations to account for strong 
electron correlation effects.!^ Plutonium based 5/ com¬ 
pounds and the isostructural samarium and ytterbium 
base d 4/ com pounds are frequently considered to be 
rpj^jJT2li3l30l3T]^ The conundrum of the TKIs arises from 
the complex interplay between many energy scales which 


involve hybridization between localized /-orbitals with 
conduction bands, strong Coulomb interactions, spin- 
orbit coupling (SOC), and hierarchy of band orderings 
driven by the crystal environment. In fact, the present 
first-principles study reveals that when we consider a 
wide range of values of U, the TKI state survives in 
SmBe, but in sharp contrast, YbBe displays a complex 
evolution of different topological phases with varying in¬ 
teraction strength. In particular, when finite U is intro¬ 
duced in the computations on YbBe, the 4/ states are 
removed from the Fermi energy (Ei?), consistent with 
recent ARPES resultJ^^^^^ the system transforms from 
being a TKI into a new d — p-type non-trivial TI in the 
intermediate range of U, and finally into a trivial band 
insulator at large value of U. 

Theoretical predictions of a TKI state in SmB^^^^^^^ 
were followed by experimental confirmations from 
transporfP^lli^, ARPE^^^H^, and scanning tunneling mi¬ 
croscopy/spectroscopy (STM/STSIn this way, the 
existence of an odd number of in-gap surface states has 
been established in SmBe. Also found is a 2D conduc¬ 
tance channel at low temperatures, giving insight into 
the puzzle of why a residual resistivity persists in SmBe 
even at very low temperatures.!^ However, theoretical in¬ 
terpretation of^ese metallic surface states still remains 
controversiaP^^^. 

Sm^m and the rare-earths YbBe and YbBiJ^ have 
also been predicted to harbor the TKI state. Early 
transport measurements indicated that YbBe is a doped 
semiconductoi!^. Recent ARPES experimenthow- 
ever find that the Yb-4/ states he ^ I eV below the 
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Ei?, suggesting a non-TKI phase, contradicting theoret¬ 
ical predictions. On the other hand, ARPES has ob¬ 
served linearly dispersing quasi-2D states without a de¬ 
tectable out-of-the-plane dispersion at the time-reversal 
invariant momenta E and X, which is consistent with the 
presence of an odd number of surface states crossing the 
Eermi level as expected in a Z 2 The temperature 

and photon energy dependence of the spectrunP^, circu¬ 
lar dichroism photoemission spectrgP^, and spin-resolved 
ARPE^^ provide further evidence to support the notion 
that YbBe is a non-trivial TL 

In this study, we examine the nature of electronic 
states in SmBe and YbBe, which share the same crystal 
structure but display very different topological phases. 
To this end, we delineate how the topological phases 
in these two materials evolve comparatively when the 
hybridization between the / and conduction levels is 
turned on or off, as well as when the strength of the 
on-site Coulomb interaction U is varied. By analyzing 
the characters of the orbitals involved, we show that 
the insulating gap in SmBe indeed stems from the in¬ 
teraction between the localized Sm-4f and itinerant Sm- 
5d bands. A parity analysis supports the presence of a 
TKI ground state. Eurthermore, the topological nature 
of SmBe is found to be essentially independent of the 
strength of electron-electron interaction. On the other 
hand, the electronic structure of YbBe is qualitatively 
similar to that of SmBe only in the small U or weak- 
coupling regime. As U is increased, the TKI phase in 
YbBe transforms into ad — p-type TI phase. In the in¬ 
termediate coupling range, the low energy bands around 
Ei? are composed of Yb-5d and B-2p orbitals, and an 
inverted gap opens up due to the SOC. With further in¬ 
crease in I/, the preceding band inversion between the 
5d and 2p orbitals is removed, and a trivial state is re¬ 
stored in YbBe but not in SmBe. We are not aware of a 
previously proposed d — p-type TI in the literature. 

The remainder of this article is organized as follows. 
In Sec. II, we discuss the crystal structure of SmBe 
and YbBe and the methodological details of our first- 
principles calculations. Section III discusses the elec¬ 
tronic structure of SmBe. Sec. IV turns to consider the 
electronic structure and the topological phase diagram of 
YbBe as a function of the strength of electron-electron in¬ 
teractions. A few concluding remarks are made in Sec. V. 


II. COMPUTATIONAL DETAILS 


Both SmBe and YbBe crystallize in the CsCl-type 
structure with the Sm or Yb ions and the Be octahe¬ 
dral clusters lying at the corners and the body centers 
of the cubic lattice, respectively, see Eig. The bulk 
Brillouin zone (BZ) is cubic, where the center of the BZ 
is the T-point, diagonal is the R-point, the center of each 
edge is the M-point, and the center of each face is the 
X-point. There are three equivalent X(M) points as a re¬ 
sult of the cubic symmetry of the crystal. We computed 
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FIG. 1. Crystal structure and the first BZ of Sm(Yb)B6. 


electro nic str uctures using the projector au gm ented wave 
met ho as implemented in the VASI®®^ package 

within the generalized gradient approximation (GGA)^ 
and GGA+lP^ schemes. Experimental lattice constants 
were used for both SmB^^and YbB^^. A 12 x 12 x 12 
Monkhorst-Pack k-point mesh was used in the computa¬ 
tions. The SOC effects are included self-consistently. 


III. ELECTRONIC STRUCTURE OF SmBe 

Fig. [^considers the evolution of electronic structure of 
SmBe as a function of U. An analysis of the three domi¬ 
nant orbital contributions to the low-lying bands reveals 
that this system is a Kondo insulator with an inverted 
band gap of 12 m eV ly ing along the F-X direction, see 
inset to Fig. The Sm-4/ orbitals clearly dom¬ 

inate the energy bands over the -0.5 eV to 0.5 eV range, 
where the bands are split into the j = 5/2 and j = 1/2 
states by the SOC (yellow dots in Fig. (a)). In the 
insulating state, the j = 5/2 states drop below the 
while the j = 7/2 states remain completely empty. The 
B-2p states lie around -1.5 eV below E^?. According to 
the ionic configuration of SmBe, Sm is a +2 ion, and thus 
the Sm-bd states are expected to be empty. In sharp con¬ 
trast, the Sm-bd states exhibit a surprisingly high degree 
of itinerancy with a bandwidth as large as 4 eV. The 
Sm-5d bands span across the Ep^ and hybridize with the 
localized Sm-4/ states near the Ep. The fractional oc¬ 
cupations of the Sm-4/ and Sm-bd orbitals results in a 
system with a mixed valence state jwhich is consistent 
with various experimental result 

Fig. shows what happens to the band structure when 
hybridization between the Sm 4/ and bd states is artifi¬ 
cially turned off. The Sm-bd states are now seen to cross 
Ep without opening an energy gap, thus recovering an 
ungapped, highly dispersive character (red dots in Fig. 
1^ (b)). The lower panel of Fig. [^further examines the 
stability of Kondo physics over a wide range of U values. 
Note that the U = 0 limit in a GGA-l-U computation 
does not imply total absence of electron correlation ef- 
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FIG. 2. (a) Band structure of SmBe from GGA calculations. 

Sizes of yellow, blue, and green dots denote weights of Sm-4f, 
Sm-5d, and B-2p orbitals in various bands, (b) Zoomed in 
view of (a) near Ef- Red dots denote metallic Sm-5d orbitals 
without hybridization with Sm-4f orbitals, (c) Low-energy 
band structure around Ef based on GGA and GGA+U (U = 
8 eV) computations. 



r 3X 3M R Gap (meV) 

Type 

SmBe (U = 

0 eV) 

- + - - 

-12 

TKI 

SmBe (U = 

8 eV) 

- + - - 

-20 

TKI 

YbBe (U = 

0 eV) 

- + - - 

-30 

TKI 

YbBe (U = 

5 eV) 

- + - - 

-3 

d-pTl 

YbBe (U = 

8 eV) 

- - - - 

+84 

Trivial 


TABLE 1. Products of parities of the occupied states at the 
eight TRIM points (T, three X, three M, and R) in various 
cases. The gap at Ef is assigned negative (positive) value for 
inverted (not inverted) band structure. 


sight into the nature of its insulating gap.^^^M^. In par¬ 
ticular, Sni-4/ orbitals are adduced to lie close to Ef, 
interacting with the highly dispersive Sm-bd orbitals as 
the tempe ratur e is lowered, consistent with theoretical 
prediction j^^ l ^^ l At the same time, the in-gap states 
are found to develop in the ARPES spectrum at low 
temperature^^^M^. Fermi surface (FS) mapping reveals 
the presence of multiple FS pockets: one pocket lies 
at F and the other two are centered on the two X- 
pointJMESl circular dichroism experiments suggest 
that the in-gap states possess chirality of the orbital an¬ 
gular momentunP^. These observations provide substan¬ 
tial evidence for the existence of topologically protected 
metallic in-gap states in SmBe- 


fects, and the mixed valence state found for U=0 can 
support Kondo physics, even though the GGA underes¬ 
timates correlation effects. Remarkably, little change in 
topology is seen between E = 0 and the case where U is as 
large as 8 eV, an increase in the hybridization gap from 
12 meV to about 20 meV notwithstanding. These re¬ 
sults demonstrate that SmBe is a Kondo insulator whose 
overall band topology is quite immune to electronic in¬ 
teractions as modeled via the on-site Coulomb energy U. 

The band structure of SmBe shows band inversion 
around the X-point in Fig. (b) with hybridization 
turned on. Since SmBe possesses inversion symmetry, the 
Z 2 topological invariants can be computed via a parity 
analysi^. The relevant products of parities of occupied 
wave functions at the time reversal invariant momentum 
(TRIM) points in the BZ are summarized in Table 1. All 
parity eigenvalues are negative except at the X-points, 
resulting in Z 2 = 1. The system is thus topologically 
nontrivial with an odd number of gapless surface states 
residing in the f-d hybridization gap. The GGA+U com¬ 
putation for U=S eV is also seen from Table I to yield 
a topological band structure like the U=0 case, albeit 
with a somewhat larger negative band gap. These re¬ 
sults agree with those of Ref. [31] in predicting that SmBe 
harbors a TKI phase. 

Recent ARPES experiments observe the presence of 
an odd number of surface states in SmBe and give in- 


IV. ELECTRONIC STRUCTURE OF YbBe 


Our computed electronic structure for YbBe is shown 
in Fig. [^for three representative values of U. At E=0, 
the band structure of YbBe resembles that of SmBe, al¬ 
though several characteristic differences should be noted. 
Firstly, unlike SmBe, where the j = 5/2 states lie around 
Ef and the j = 1 {2 states are completely empty (Fig. 

(a)), the 2 — 5/2 states of Yb-4/ orbitals are fully oc¬ 
cupied in YbBe (-1.5 eV below Ef), while the j = 7/2 
states lie near Ep (Fig- (a)). Secondly, the insulat¬ 
ing band gap is larger in YbBe with a value of 30 meV 
along the X-M direction, see inset to Fig. (d). On 
the other hand, similar to the case of SmBe, the narrow 
energy gap in YbBe stems from the hybridization be¬ 
tween Yb-4/ and Yb-5d orbitals, and the parity analysis 
of Table I reveals that the system is a TKI for small val¬ 
ues of U (Fig. (d)).l^ Experimental data , how ever, do 
not support these theoretical predict ions Recent 

ARPES experiment^^^Mm show that Yb-4/ states are lo¬ 
cated at ^ -1 eV and -2.3 eV below Ef, indicating that 
YbBe does not belong to the Kondo insulator paradigm. 

As we turn on electronic interactions via GGA+U cal¬ 
culations, we see substantial changes in the band topol¬ 
ogy of YbBe in the two right columns of Fig. Yb- 
4/ bands drop below Ef as we increase the value of E, 
with two narrow Yb-4/ bands located at -0.5 eV (-1.4 
eV) and -1.8 eV (-2.7 eV) binding energy for U = 5 
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FIG. 3. Band structure of YbBe for (a) = 0, (b) = 5 eV, 

and (c) = 8 eV. Black solid lines and light orange circles 

in (a) give GGA+SOC and GGA+SOG+GW results, respec¬ 
tively. (d-f) Zoom in views of band structures around the Ep 
corresponding to (a)-(c). Sizes of yellow, blue, and green dots 
give weights of Yb-4/, Yb-5d, and B-2p orbitals in the band 
structures, (g-i) Schematic band structures corresponding to 
the results of panels (a)-(c). 


eV (U = 8 eV). As Yb-4f bands move to higher binding 
energies, their hybridization with conduction bands be¬ 
comes smaller. The B-2p orbitals, on the other hand, are 
pushed toward Ep with increasing U and interact with 
Yb-5d states (Fig. (e)). The low energy band struc¬ 
ture clearly shows a band inversion between the Yh-bd 
and B-2p bands. Combined with the inversion of orbital 
texture, the parity analysis given in Table I indicates 
that YbBe is a non-trivial TI for U values around 5 eV. 
YbBe thus presents a no vel TI state involving inversion 
of d-p type orbit als.l^^J^ We have recalculated the band 
structure using a GGA+SOC+GW scheme, and found 
results similar to those obtained via GGA+SOC calcula¬ 
tions, indicating that d-p type band inversion cannot be 
obtained via the GW correction alone, see Fig. [^(a). 

Our preceding scenario is in substantial accord with 
ARPES data on YbBe, which exhibit nearly linearly dis¬ 
persing quasi-2D states without a significant out-of-the- 
plane momentum (k^) dispersioiPSlM]^ SmBe, the 

FS of YbBe consists of an oval-shaped contour and a 
nearly circular contour enclosing X- and F-points. The 
temperature evolution of the states does not reveal a 
Kondo-like behavioJ^. Circular dichroisnP^ and spin- 
resolved ARPE^ spectra show that the surface states 
are spin-polarized and display spin-momentum locking. 


These experimental results support YbBe being a TI, 
but not a TKI. 

With further increase of t/, not only the Yb-4/ bands 
move to higher binding energies, but the band gap be¬ 
tween Yh-bd and B-2p states closes and reopens, and the 
band inversion at the X-points is removed, see Fig. 

Signs of parity products at all TRIM points are the same 
for U>6.5 eV, and the system thus assumes a trivial band 
insulator ground state at large U. 
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FIG. 4. (a) Weights of Yb-4f (red dots) and B-2p (blue dots) 

orbitals at the bottom of the conduction band at the X-point 
as a function of U in YbBe. Lines are drawn through the 
computed points to guide the eye. (b) Gorresponding direct 
band gaps at the X-point. 

Fig-i gives insight into the key reasons driving the 
existence of two distinct nontrivial topological phases in 
YbBe. Fig.j^a) analyzes the contributions of Yb-4/ and 
B-2p orbitals as a function of U at the bottom of the 
conduction band at the X-point where the band inversion 
occurs. With increasing U, weight of the Yb-4f orbitals 
decreases, while that of B-2p orbitals increases, indicat¬ 
ing that Kondo physics weakens near the Ep. The or¬ 
dering between the Yb-4f and B-2p levels reverses around 
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U ^4 eV, and the B-2p orbitals start to dominate near 
Ef> This characteristic evolution can also be seen in the 
band gap in Fig. [^b). The band gap at Ep initially in¬ 
creases weakly with increasing U as long as the gap is 
defined between the Yb-derived / and d bands. Around 
U ^3 eV, the insulating gap switches over to become a 
gap between the Yb-5d and B-2p states, and decreases 
sharply for higher U values. The system thus under¬ 
goes a smooth transition from being a TKI to a d — p 
type TI without closing the gap. This novel topologi¬ 
cal phase transition is consistent with adiabatic continu¬ 
ity arguments, which dictate that the topology of two 
band structures is the same, if one can go from one to 
the other through a series of adiabatic transformations 
of the Hamiltonian without closing the gap anywhere in 
the BZ. YbBe loses its topological character as the band 
gap disappears around U ^6.5 eV, and reopens to yield 
the trivial band insulator at higher values of U. In this 
way, modulation of Hubbard U can vary the topological 
state and lead to a rich phase diagram.^ Tuning of the 
Coulomb interaction in a given material is of course not 
straightforward to realize experimentally. It could per¬ 
haps be accomplished, however, via doping with Sm, Pu 
or other elements with the same crystal structure, which 
may change the value of effective U through changes in 
Coulomb screening, and provide a pathway for verifying 
our predictions. 

V. CONCLUSIONS 

We have systematically investigated the topological 
phase diagram of mixed valence compounds SmBe and 
YbBe by varying the strength of on-site Coulomb inter¬ 
action U within the framework of the GGA-l-U scheme. 
The hybridization gap in SmBe between the localized 
Sm-4f bands and the conduction bands remains fairly 
unchanged as U is varied up to values as large as 8 eV. 
In sharp contrast, YbBe exhibits a much more complex 
behavior. As U is increased, Yb-4/ states in YbBe grad¬ 
ually move to lower energies, while the weakly correlated 
B-2p bands move upward, and these two bands eventu¬ 


ally get interchanged without closing the insulating band 
gap around the Ep- In this way, the TKI phase assumed 
by YbBe for small (weak coupling) values of U trans¬ 
forms into a new d-p type of topological phase for values 
of U in the intermediate coupling regime. The topologi¬ 
cal phase in YbBe reverts to a non-topological state for 
large values of U, but this does not happen in SmBe. 
Our results concerning the topological behavior of pris¬ 
tine SmBe and YbBe are in reasonable accord with the 
available ARPES results on these compounds, and give 
new insight into the comparative evolution of their band 
topologies with increasing interaction strength U. Our 
theoretical predictions concerning the novel topological 
phases of YbBe might be amenable to experimental test¬ 
ing by substituting Yb by Sm, Pu or other similar el¬ 
ements to drive changes in the value of effective U via 
changes in Coulomb screening effects in the material. 
Due to the involvement of d electrons in its topologi¬ 
cal phase, YbBe would be prone to quantum fluctuations 
and instabilities, offering unique opportunities for inves¬ 
tigating the robustness of topological protections in the 
presence of such instabilities. 
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